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ABSTRACT 

In this paper, the question of the motion of a projectile 
thrown at an angle to the horizon is learned. With zero 
air drag force, the analytic resolution is well 
understood. The course of the projectile is a parabola. 
In situations of practical interest, such as jaculatory a 
ball with the event of the strike of the medium the 
square resistance law is generally used. In that casing 
the question likely does not have an true analytic 
resolution and therefore in most expert publications it 
is solution numerically. Analytic advance to the 
resolution of the question are not sufficiently 
professional. Meanwhile, resolve solutions are very 
handy, for a honest preparedness to win problems and 
are expressly costly for a qualitative analysis. That is 
why the representation of the missile summon with a 
uncombed rough resolve formula under the quadratic 
air resistance instant big methodological interest. 
Lately these formulas have been procure. These 
formulas bestow us to obtain a finished divisive 
narration of the problem. This delineation perfect 
divisive formulas for shape the basic eight parameters 
of missile movement. Analytical formulas have been 
infer for the six fundamental official dependences of 
the statement including the trajectory equality in 
Cartesian coordinates. Also this description contains 
the judgment of the optimal throwing angle and 
maximum frequent of the motion. In the destitution of 
publicity resistance, all these relations transfer into 
well-understood formulas of the h 5 q)othesis, 
speculation of the parabolic motion of the projectile. 
The proposed analytical discharge oppose from other 
solutions by clearness of formulas, comfortableness of 
use and high fidelity (relation hallucination is near 1-2 
%). The motion of a baseball is instant as an example. 


The converse formulas make it possible to carry out 
an resolvent examination of the motion of a projectile 
in a medium with resistance in the highway it is done 
in the casing of no drag. 

Keywords: spherical object, quadratic drag force, 
analytical formula, relative error 

Introduction 

The problem of the motion of a detail mass 
(projectile) thrown at an angle to the circle in midair 
has a g history. It is one of the great humanistic 
problems that lead the dynamics as knowledge. The 
amount of works devoted to this employment is vast. 
It delineatea elemental of many prefatory courses of 
physics. With no trail, the analytic resolution of this 
question is well known. The course of the stage 
assemble is a parabola. In this notes, an analytic 
appropinquate is necessity for the scrutation of the 
projectile motion in a medium with tetragonal 
resistance. The opposed of the immediate product is 
to give simple formulas for the sense of the projectile 
trajectories under the course with quadratic 
intelligence resistance. These formulas are 
advantageous even for first-year undergraduates. The 
question of the motion of a projectile in midair has 
stirred up interest of authors (Cohen et al., 2014; 
Kantrowitz & Neumann, 2013; Borghi, 2013). For the 
construction of the separative solutions various 
methods are application - both the traditive 
approximate (Benacka, 2010; Vial, 2007; Parker, 
1977; Erlichson, 1983; Tan, Frick & Castillo, 1987)], 
and the modem methods (Yabushita, Yamashita & 
Tsuboi, 2007). All proposed approximate analytical 
solutions are rather complicated and uncomfortable 
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for educational aim. This is why the description of the 
projectile course by means of a single approximate 
anal54ical formula under the square mien resistance is 
of powerful methodological and instructive 
significance. In (Chudinov, 2002; 2004; 2013) 
comparatively harmless approximate analytical 
formulas have been possess to ponder the motion of 
the projectile in a medium with a quadratic comfit 
violence. In this article, these formulas are interest to 
solve the humanistic problem of maximizing the 
projectile distance. From now on, the word “appoint 
mass” contemptible the center of mass of a smooth 
rounded object of limited circle r and trial-sectional 
region S = 7ir2. The provision of applicability of the 
quadratic resistance law are esteem to be fulfilled, i.e. 
Re 5 molds scalar Re sleep within 1x103 < Re < 2x105. 
These values correspond to the missile motion 
velocity, fabling in the range between 0.25 m/s and 53 
m/s. 

Equations of Point Mass Motion and Analytical 
Formulas for Basic Parameters 

Suppose that the force of gravity concern the appoint 
mass together with the force of air resistance R (see 
Fig. 1). Air resistance force is corresponding to the 
square of the celerity of the point body and is 
addressed opposite the fleetness vector. For the 
convenience of further calculations, the starch cascade 
will be written as 2 R mgkV= . Here m is the mass of 


the projectile, g is the acceleration due to gravity, k is 
the proportionality agent. Vector equation of the guide 
of the point mass has the beauty 


mw = mg + R, 


where w - acceleration vector of the point mass. 
Differential equations of the motion, commonly used 
in ballistics, are as follows (Okunev, 1943) 


i!^=-gsine-gkV^ = il = Vcose 

dt ^ ^ ’dt V ’ dt 

= v sine (1) 

’dt ^ ^ 


Here V is the velocity of the appoint sum, 0 is the 
angle between the tangent to the course of the 
moment bulk and the flat, x, y are the Cartesian 
coordinates of the instant bulk, k is 


^ _ PaCdS 

2 mg 


1 

-77 = conts, 


Pais the air density, cd is the drag factor for a sphere, 
S is the misfortune-section scope of the goal, and 
V_term is the termination celerity. The first two 
equations of the system (1) personate the projections 
of the vector equation of guide on the tangent and 
main analogical to the course; the other two are 
kinematic relations connecting the projections of the 
fleetness vector point bulk on the shaft x, y with 
derivatives of the coordinates. 


Figure 1: Basic Motion Parameters 
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The well-known solution of system (1) consists of an 
explicit analytical dependence of the velocity on the 
slope angle of the trajectory and three quadratures 


v(e) = 


VO cose 


cosO Jl+kVQCOs^0o( f (^o)“/(^) 


fm = 


sin0 

cos20 


+ in 



( 2 ) 


T=to--fa -^dd.x 

g ■’Oo COS^ 

y=yo-^Qv^tgddd 



( 3 ) 


Here VO and 0 are the commencing importance of the 
fleetness and the slope of the trajectory respectively, 


to is the incipient importance of the period, xO, yO are 
the drop cap values of the coordinates of the instant 
bulk (most often approve t_0=x_o=y_0=0). The 
derivative of the formulas (2) is shown in the well- 
known monograph (Timoshenko & Young, 1948). 
The integrals on the right-hand sides of formulas (3) 
cannot be expressed in limit of simple performance. 
Hence, to shape the variables t, x and y we must either 
complete system(l) numerically or appraise the 
precise integrals (3). Comparatively single 
approximate analytical formulas for the principal 
parameters of tendency of the projectile are deduce in 
(Chudinov, 2002, 2004). The four parameters address 
to the top of the trajectory, four - to the stage of drip. 
We will give a consummate concise of the formulas 
for the maximum height of ascent of the point mass 
H, motion time T, the fleetness at the trajectory tip Va 
= V(0), stampede range L, the repetition of ascent ta , 
the absciss of the trajectory pinnacle xa, impact angle 
with esteem to the horizontal 01 and the final speed 
VI (see Fig. 1). These formulas are sum in the right 
column of Table 1. In the port column of this Table 1 
homogenous formulas of the allegorical theory are 
immediate for illustration. 


Analytical Formulas for the Main Parameters 


No drag 
(R = 0) 

Quadratic drag force ( R = mgkV 2 ) 

^_visiv}eo 

V^sin^Bo 

5(2+fcvf‘""o) 

OsinGQ _ 2 \2H 

9 9 


Va — l^OCos^O 

FnCOSSn 

V - u U 

+ kPo^ l^sin^o + cos^Gq. in tan 

l=^V^sin2eo = VJ 

L=FaT 

VoSindo T 

g 2 

T - kHVa 
ta= 2 

II 

f (0i) 


With zero drag (k = 0), these formulas go over into 
the wary formulas of the point mass parabolic motion 
theory. All movement characteristics relate by these 
formulas are functions of opening requisite of 
throwing VO, 00. Proposed formulas have a bounded 
region of application. We present the symbolology 2 


Op kV_0^2 . The dimensionless parameter p has the 
profession physical object - it is the ratio of air 
resistance to the moment of the projectile at the 
opening of the motion. The principal characteristics of 
the motion H, T, Va , L, xa have precision to within 2 
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- 3% for values of the pierce angle, for commencing 
speed and for the parameter p from ranges 

0°< 00 < 70°, 0 0 50 V □ □ m/s , 0 < p < 1.5. 


calculating the main parameters. Now it is possible to 
construct the trajectory in the entire range of launch 
angles and at values of the initial velocity and the 
parameter p 


For a baseball the typical values of the drag force 
coefficient k are near 0.0005 ^ 0.0006 s2/m2, 
maximal commencing speed is circularly 50 m/s 
(Cohen et al., 2014; Tan et al., 1987). Therefore the 
converse formulas are competent for the qualitative 
and quantitative narrative of the motion of the 
baseball and other resemblant goal. These formulas, 
in devote, occasion it likely to obtain a 
uncompounded divisive formula for the principal 
functional relationship of the proposition y(x) 
(Chudinov, 2002). In the want of air resistance, the 
trajectory of a point mass is a parabola. The equation 
of the trajectory can be scriptory in two constitution. 
It can be scriptory in bound of the opening station of 
jaculatory Vo,^o (first figure). It can also be written in 
terms of the motion parameters H, L, xa (second 
form) 


Y(x)=x.tan0o ~ 


gx^ 

2v^cos^^o 


Hx(L-X) 



( 4 ) 


The course is symmetric with respect to the greatest. 
When the item mass is under a drag force, the course 
becomes asymmetric. The top of the trajectory is 
divide towards the point of incident. In title, a vertical 
asymptote appearance near the trajectory. Taking 
these circumstances into narrative, the performance 
y(x) may be arrange using parameters H, L, xa as 
(Chudinov, 2002) 


Y(x) 


. HxiL-X) 
'xl+(L-2Xa)x 


( 5 ) 


The constructed dependence y(x) provides the shift of 
the apex of the trajectory to the right and has a 
vertical asymptote. In the case of no drag L = 2xa 
and formula (5) goes over to formula (4). We note the 
remarkable property of formula (5). We insert the 
exact values of the parameters L, H, obtained by 
numerical integration of system (1), into formula (5). 
Then the numerical trajectory and the analytical 
trajectory constructed by means of formula (5) are 
identical. This means that formula (5) approximates 
absolutely precisely the projectile's trajectory which is 
numerically constructed using system (1) at any 
values of the initial conditions • The some 

transformation of the proposed formulas (Chudinov, 
2013) makes it possible to improve the accuracy of 


0°< 6>o< 90°, 0< Vo ^ 80 m/s , 0 < p < 4 . 

As an represent of the use of the speak formulas from 
the Table 1 and of formula (5) we calculated the 
movement of a baseball with the following initial 
station 

Vo = 40 m/s ,00 = 45°, k = 0.000625 s^/m^ g = 
9.81 m/s^ 


Comparison of Numerical and Analytical Calculations 


Parameter 

Analytical 

value 

Numerical 

value 

Error 

(%) 

H(m) 

30.1 

29.8 

1.1 

T(s) 

4.96 

4.91 

1.0 

Va (m/s) 

19.3 

19.3 

0 

L (m) 

95.7 

96.1 

-0.4 

ta (s) 

2.30 

2.31 

-0.4 

xa (m) 

53.7 

53.0 

1.3 

01 (degree) 

-58.6° 

-57.3° 

2.2 


The terminate of calculations are recorded in Table 2 
and are shown in Figure 2. The help atlantes restrain 
the importance calculated by deliver formulas from 
the Table 1. The third atlantes show the importance of 
parameters obtained by numeral integration of the 
motion equations (1) by the arrangement RungeKutta 
process. The deviations from the correct values of 
parameters are shown in the fourth column of the 
Table 2. Tabulated data show that the importance of 
basic parameters of the projectile motion (flight range 
L, tendency time T, prominence H) adapted by 
divisive formulas contend from the exact values no 
more than 1%. An round trajectory (5) is compile with 
the help of parameters H, L, xa. It is shown in Figure 
2 (dotted line). The indistinct important rope in Figure 
2 is possess by numerical integration of system (I) 
with the befriend of the 4-th order Runge-Kutta 
method. As can be accomplished from Figure 2, the 
divisive solution (formula (5)) and a numeral 
resolution are almost the same. The frustrate line in 
Figure 2 is erect in the destitution of mien resistance 
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The Graph of the Trajectory y= y(x) 


Formula (5) 


Y, m 



One of the most necessary aspects of the projectile 
movement problem is determination of an optimal 
viewpoint of jaculatory of a point number which 
provides the maximum range. Let us make usefulness 
of formulas from the Table 1 and (5) to resolve the 
problem of optimization of motion parameters 


Determination of an Optimum Throwing Angle and 
Maximum Range 

Let the point of incidence A be on a horizontal 
straight line defined by the equation const y= 
yi=const 


Statement of the Optimization Problem 



To solve the problem, let us companion use of the 
appoint mass trajectory equality (5). In this equality, 
the tendency characteristics H, L, Xa are duty of initial 
station of throwing, Vo^o- They are explain by 
formulas from the Table 1. Let us change yi in the 
sinistral-hand part of equality (5) and resolve it for 
variablex. Thus, we must get the fleeing frequent 
formula for the casing under contemplation 


SUMMARY 

The proposed approach based on the usefulness of 
analytic formulas mate it possibility to disintricate 
way a qualitative analysis of the motion of a projectile 
with the air draw taken into reckoning. All basic 
parameters of guide and various problems of 
optimization are relate by unadorned analytical 
formulas confine rudimental functions. Moreover, 
numeral values of the sought variables are resolute 
with an comfortable propriety. It can be accomplish 
even on a authoritative calculator. Lately some creator 
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(Warburton & Wang, 2004; Stewart, 2005; Hu et al., 
2012) have used the Lambert W service to study the 
projectile motion with resistance. But this relatively 
“new” function is not profitable on a calculator. 
Special algorithms are required to compute this 
service. Thus, discourse formulas make it possibility 
to study projectile motion with tetragonal drag force 
even for first-year undergraduates. In end, we hope 
that efforts to prevail an analytical release to this 
question will be extended and will achieve new exact 
and effective solutions. 

CONCLUSION 

A projectile that does not contain an cordite charge or 
any other good of instruct is condition a kinetic 
projectile, motory action weapon, seepy Life warhead, 
dynamic warhead or motory penetrator. Typical 
kinetic energy weapons are blunt projectiles such as 
cliff and round shots, pointed ones such as directed 
edge, and little pointed ones such as ball. Among 
projectiles that do not hold explosives are those 
launched from railguns, coilguns, and quantity 
drivers, as well as active force penetrators. All of 
these weapons work by obtain a lofty muzzle velocity, 
or commencing velocity, collectively up to 
(hypervelocity), and collide with their aim, translate 
their kinetic energy into mischievous shock waves 
and heat. Other types of animated arm are accelerated 
over era by a rocket engine, or by gravity. In either 
plight, it is the active strength of the projectile that 
consume its target. 
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